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Abstract 

This paper examines the properties of a regularization of Burgers 
equation in one and multiple dimensions using a filtered convective ve- 
locity, which we have dubbed as convectively filtered Burgers (CFB) 
equation. A physical motivation behind the filtering technique is pre- 
sented. An existence and uniqueness theorem for multiple dimensions 
and a general class of filters is proven. Multiple invariants of motion 
are found for the CFB equation and are compared with those found 
in viscous and inviscid Burgers equation. Traveling wave solutions are 
found for a general class of filters and are shown to converge to weak 
solutions of inviscid Burgers equation with the correct wave speed. 
Accurate numerical simulations are conducted in ID and 2D cases 
where the shock behavior, shock thickness, and kinetic energy decay 
are examined. Energy spectrum are also examined and are shown to 
be related to the smoothness of the solutions. 
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1 Introduction 



The Euler and Navier- Stokes equations are well known as the fundamen- 
tal laws governing fluid dynamics, yet even after 200 years they continue 
to present theoretical and computational challenges. The nonlinear terms 
inherent in the equations give rise to small scale structures, in the form of 
turbulence and shocks, which have proven to be the bane of computational 
simulations. With the proper modeling of these small scales, we hope it is 
possible to address both turbulence and shocks with one encompassing tech- 
nique. Currently, the Lagrangian averaging approach is making strides in 
handling turbulent flows PElElHlElEllSElElliniin]. That work moti- 
vated the work presented in this paper. The Lagrangian averaging approach 
results in a filtered convective velocity in the nonlinear term. This paper also 
uses a filtered convective velocity in the nonlinear term of Burgers equation, 
with the intention of discovering if this technique is a reasonable means of 
capturing shock formation. 

Burgers equation has been a useful testing grounds for fluid dynamics 
for many years due to the fact that it shares the same nonlinear convective 
term as Euler and Navier-Stokes equations. Indeed, Burgers equation has 
been the focus of much work, both numerically and analytically [TJl H31 HH 
[TBI ITS] ITT] [T5| \W[ [2D] . Originally Burgers equation was introduced and used 
as a simplistic model for one dimensional turbulence (J2J 02] • The original 
equation, 

u t + uu x = UU XX , (1) 

which shall be referred to as viscous Burgers equation, includes a dissipative 
viscous term. By removing the viscous term, 

Ut + uu x = 0, (2) 

inviscid Burgers equation is obtained. Much like the Euler equations, inviscid 
Burgers equation can be expressed as a conservation law. Both viscous and 
inviscid Burgers equation are easily extended into multiple dimensions giving 

u t + u ■ Vu = v A u (3) 

and 

U( + u-Vu = 0. (4) 

It is well established that inviscid Burgers equation forms discontinuities 
in finite time, determined by initial conditions [T5| [16]. To deal with these 
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discontinuities weak solutions are introduced. However, when weak solutions 
are introduced, solutions are no longer necessarily unique p21 EI] . In order 
to choose the physically relevant solution, an entropy condition is applied, 
which one and only one weak solution satisfies. This physically relevant 
solution is referred to as the entropy solution. Lax, Oleinik, and Kruzkov 
have examined the entropy condition for conservation laws and expressed it 
using different techniques [TU EH E2] • Each of their entropy conditions can 
be used in different classes of conservation laws, but can all be applied to 
inviscid Burgers equation with equivalent results [23J. 

With the dissipative term, viscous Burgers equation does not form discon- 
tinuities [T51 US]. Indeed it smooths discontinuities found in initial conditions, 
and has a unique infinitely differentiate solution for all time [El EI] . More- 
over, the limit of the solutions as v — > converge strongly to the entropy 
solution of inviscid Burgers equation [161 EU 122] . 

Viscous Burgers equation is not the only regularization of inviscid Burgers 
equation that converges to the entropy solution. Many regularizations have 
been proposed, typically with the addition of a dissipative term. Among them 
are regularizations with hyperviscosity, filtered viscosity, and a combination 
of viscosity and dispersion [T9l 1251 1251 [271 128] . 

Another famous regularization is the KDV equation which uses a disper- 
sive term, 

u t + uu x = -eu xxx . (5) 

This regularizes inviscid Burgers equation in the sense that solutions are now 
continuous, however, many oscillations form as e — > 0, requiring a weak limit 
for convergence [251 ED] • This limit is not the entropy solution, nor even a 
weak solution of inviscid Burgers equation [31] . 

The incompressible isotropic LANS-a equations are given by [10] 



— + u-Vu + u J -V% = -VP+— -Au (6) 
ot Re 

where u is defined as 

u = u — a 2 A u. (7) 

As mentioned earlier this equation has had success capturing some of the 
behavior for a class of turbulent flows. It can be seen in Equation (|6| that 
LANS-a employs a filtered convective velocity. What this paper proposes 
is that it is this filtered convective velocity in the nonlinear term that is 
successful in modeling small scale behavior. If this is indeed the case, then 
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such a term could also be used to model shocks alongside turbulence. This 
paper begins the examination of using such a nonlinear term in shock regu- 
larization. As was stated, Burgers equation has been extensively researched 
and is well understood. It forms shocks readily and has been chosen as a 
testing ground for using a filtered convective velocity for shock regulariza- 
tion. This technique is not intended to be an analytically or numerically 
superior method of regularizing shocks, but more as a proof of concept with 
the intention of extending this method into areas involving both shocks and 
turbulence. 

This paper considers what will be referred to as convectively filtered Burg- 
ers equation (CFB), where the convective velocity in inviscid Burgers equa- 
tion is replaced with a filtered (averaged) velocity. The filtering is done by 
convoluting the velocity u with the kernel g (properties of which will be 
discussed in section [3]) resulting in 

Ut + uu x = (8) 
u = g * u. (9) 

If an inverse for the filtering exists, the equation, using only the filtered 
velocity, can be rewritten as 



u t + uu x = g 1 (uu x ) - ug 1 (u x ), (10) 

where g -1 represents the inverse of the low pass filter. Similar to the Large 
Eddy Simulation (LES) or the Lagrangian Averaged Navier-Stokes-a (LANS- 
a) equations where the effect of small scales (high wavenumbers) on the large 
scales are modeled by a subgrid scale stress tensor [321 EH] • To this end the 



right hand side of Equation (10) can be considered a convective subgrid 
scale stress for regularizing a shock. 

This being said, it is the intention that the filtered velocity, u, be consid- 
ered the physically relevant quantity, that captures the proper macroscopic, 
or low wavemode, behavior. In Sections [7 10, u from CFB will be compared 



to the velocity from viscous Burgers equation. The unfiltered velocity, u is 
intended to capture both the low and high wavemode activity. In Section [5] 
most of the invariants of motion found for CFB involve the unfiltered velocity. 

CFB, like Burgers equations, is easily extended into multiple dimensions, 
with equation 



4 



u t + u-Vu = (11) 
u = g * u. (12) 

Special attention is paid to the Helmholtz filter, which gains notice for 
several reasons. It is a common filter, employed in the Leray-a model of 
turbulence [HEJE], Lagrangian Averaged Navier-Stokes (LANS-a) [UE1IM1 
El E| and in the Lagrangian averaged Euler-a 0, [TO]. It has been shown 
that one dimensional CFB is Hamiltonian when the Helmholtz filter is used 
[351 136] . Also, the inverse of the Helmholtz filter is known for vanishing or 
periodic boundary conditions, so it can be expressed through convolution or 
using differential notation. In ID, the following are equivalent 



u = ^e~' '' * u (13) 
u = u- u xx . (14) 

When moving to higher dimensions the convolution kernel changes and 
the double derivative becomes a Laplacian. Using the Helmholtz filter, Equa- 



tion (10) can be simplified to 

u t + uu x = —3a 2 (I — a 2 dx 2 )~ (u x u xx ). (15) 

The CFB equation is not the first look into the use of filtered convective 
velocities. Jean Leray proposed using a filtered convective velocity in the 
Navier-Stokes equations as early as 1934 [37]. This was done in an attempt 
at proving that regular solutions to the Navier-Stokes equations exist and to 
examine properties of those solutions. This has led to several investigations 
into the Leray-a model of turbulence [TJ [38j [39] primarily using the Helmholtz 
filter. In the projected models for extending our technique into the Euler 
equations, the filtering occurs in more than just the convective velocity. At 
that point the research loses much of the similarity with Leray's work. This 
extension into Euler equations will be the topic of following papers. 

The CFB equation itself has been previously studied. The family of 
equations 

Ut + uu x + buu x = 0, 
u = g * u 
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where b was a free parameter, has been examined [31)1 HP] , where g was almost 
exclusively the Helmholtz filter. The b = case reduces to the ID version 
of CFB. Holm and Staley [35J established the invariant quantities, J u and 
II^Hoo. They examined peakon solutions and cliff solutions both analytically 
and numerically. A Lagrangian structure was shown for the b ^ case and 
a Hamiltonian structure proposed for b = by Degasperis et al. [ID]. Bhat 
and Fetecau [36] also examined one dimensional CFB with the Helmholtz 
filter. They more fully established the Hamiltonian structure proposed by 
Degasperis et al. They also proved existence and uniqueness of a solution, 
and proved that the solutions converge to a weak solution of inviscid Burgers 
equation. 

This paper extends the results into the use of more general filters and 
looks into ID and multiple dimensions. It also examines many of the physi- 
cally relevant characteristics of CFB. Some of the results in this paper have 
been previously presented in conference papers [H HU, H2] . The next section 
elaborates on the motivation for the averaging (filtering) technique. Section[3 
establishes the characteristics desired in the filter. Existence and uniqueness 
theorems are given in section [4] for CFB in multiple dimensions and a general 
class of filters. Section [5] examines the various invariants of motion for CFB 
and compares them to those found in viscous and inviscid Burgers equation. 
In section [6j a traveling wave solution for general filters is presented and an- 
alyzed. Sections [7 10 deal with numerical simulations and results regarding 
shock behavior, spectral energy, and energy norms. All is then followed with 
concluding remarks. 



2 Shock Formation 

We begin by examining the mechanics behind shock formation. The non- 
linear convective term u ■ Vu generates high wave modes, by transferring 
energy into smaller scales as time progresses. This term can be found in 
the Euler and Navier-Stokes equations, where it is responsible for tilting and 
stretching of vortical structures [IB]. It is also found in Burgers equation. If 
this cascade of energy into the smaller scales is left unchecked, discontinuities 
will form. To prevent this, viscosity is added. Dependent on the Reynolds 
number, there exists a length scale where viscosity begins to dominate the 
energy cascade, transferring kinetic energy into thermal energy through dis- 
sipation. In the Navier-Stokes equation, the energy cascade has a slope of -| 
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until this length scale, known as the Kolmogorov scale, illustrated in figure 
IJi. Viscous Burgers equation has an energy cascade slope of -2 until viscos- 
ity begins to dominate; seen in figure [TJj [1U EH] . The Euler equations and 
inviscid Burgers equation are the respective limits as viscosity goes to zero. 
By filtering the convective velocity, the energy cascade is affected through 
the nonlinear term. A low pass filtered velocity will have less energy in the 
high wave modes, thus the nonlinear term u ■ Vu will generate higher wave 
modes at a reduced rate. By reducing the energy cascade, discontinuities are 
prevented from forming. 




Figure 1: Energy cascade sketches for Navier-Stokes/Euler equations and 
viscous/inviscid Burgers equation, (a) Energy cascades from high wave- 
lengths to lower wavelength at a predicted rate of -| for Navier-Stokes/Euler 
equations. Kinetic energy drops drastically upon reaching a certain wave- 
length. That wavelength 77 is referred to as the Kolmogorov scale. Here 
Rei < Re 2 < Re 3 . (b) For Burgers equation the energy cascade has a slope 
of -2 until viscosity begins to exert its influence. Here u 3 < < v\. 



3 Filters 



In Equation ( 12 ) , the filtered velocity is obtained by convoluting the unfil- 
tered velocity u with a low pass filter g. This section discusses the require- 
ments and properties of the low pass filter. Using physical and analytical 
arguments the class of relevant filters is outlined. 



3.1 Filtering as a Weighted Average 



From a physical perspective, Equations (11) and (12) describe a fluid 
where the convective velocity of a particle is the weighted average of the 
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velocities of the particles around it. In such an averaging, there are several 
guidelines that seem intuitively reasonable. The averaging should give no 
particle a negative weight, should give no directional preference, and should 
give more weight to particles that are physically closer. These guidelines are 
easily expressed as mathematical properties of the filter. Thus g should be 
nonnegative, decreasing, and radially symmetric. Furthermore, as a weighted 
average the integral of g over the domain should be precisely one. 



3.2 Characteristic Wavelength Parameter 

Each such filter g is equipped with a parameter a, such that a is the char- 
acteristic wavelength of that filter. This parameter is introduced by scaling 
the filter as such: 

f - y-). (16) 
This scaling is also realized in the Fourier domain by noting that 

P(k)=£(ak). (17) 

Thus as a becomes smaller, the wavelength where the filter exerts influence 
also become smaller. 

With this scaling, the filter remains normalized, nonnegative, decreasing, 
and isotropic. Furthermore as a — > 0, g a approaches the Dirac delta func- 



tion, and Equation (11) formally limits to inviscid Burgers equation. This 
parameter is crucial in resolving features in u. When convoluting g a with u, 
features in u that have length less than a will be smoothed out. It is also 
seen in Section [6] that a controls shock thickness. 



3.3 Fourier Decay Properties 



In the following section, Theorem 4.1 states that system (21) , (22) 
and (23) has a continously different iable solution, if g(x) e W 1 ' 1 



alternatively 



and 



\g\W < °°> 



d 



g\\ L i < oo. 



l ) or 
(18) 

(19) 
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Properties 



Mathematical Expression 



Normalized 
Nonnegative 
Decreasing 
Symmetric 
Fourier Decay 



1 9=1 
<?(*) > 0, Vx 

|xi| > |x 2 | flf(Xi) > #(x 2 ) 
|xi| = |x 2 | flf(xi) = flf(x 2 ) 
limik^oo |k|g(k) = 



Table 1: This table succinctly lists the requirements of the low pass filters 



employed in Equation (12) 



Condition (|19|) has implication in the Fourier domain. The Riemann- 

ikj§ G Cq. Thus 



Lebesgue Lemma demands that if -$-g G L l then 
not only must g & C , but the condition, 



|k| 



lim |k|£(k) = 0, 



(20) 



must be satisfied in order to meet (19). Note that ( |20| ) is a necessary, not 
sufficient, condition to meet Theorem 4.1 



4 Existence and Uniqueness Theorem 



It has been previously proven [36J that the initial value problem (21) , (22) 
, and (23) using the Helmholtz filter (13) has a continuously differentiable 



solution for continuously differentiable initial conditions. Taking inspiration 
from that work, the following theorem generalizes the existence and unique- 
ness result into multiple dimensions and a variety of filters. 



Theorem 4.1 Let #(x) G W 1 ' 1 (R n ) 
a unique global solution u(x, t) G C 1 



and Uq(x) G C 1 ) 



then there exists 



to the initial value problem (21) , 



(22) , and (25) . 



u t + u • Vu = 
u = g * u 
u(x,0) = u (x) 



(21) 
(22) 
(23) 
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Proof Begin by shifting perspective into the material view. Associate a map 
</>(£,£) : M. n — > W 1 as the map from a particle's original position, (£), to the 
particle's position at time t, (x). Naturally this dictates that 0(^,0) = £. 
Associate this material map with the velocity u by 



d_ 

dt 



4>(&t) = u(0(£,t),t). 



It can then be seen that (21) simply becomes 



d_ 
dt 



u(0(£,t),t) = 0, 



and thus, u(0(£,t),t) = u(0(£,O),O) = u (£), which implies 



u 



u 



(24) 



(25) 



(26) 



If it is assumed that 0(£,i) has a continuously differentiate inverse 
_1 (x,t), then (21) , (22) , and (23) will have the continuously dif- 
ferentiate solution 

u(x,t) = u o (0- 1 (x,t)). (27) 

A sufficient condition for such an to uniquely exist is the Jacobian of 
4> to be non-zero for all positions and time. Thus if J(<p($,)) ^ 0, V£,t, then 
(@ , (g , and (g is uniquely solved by ^ . 

Using a result from Aris [16] , the time derivative of the Jacobian is found 
to be 



d_ 
dt 



J = V • u J, 



(28) 



(29) 



which is a different iable equation with solution 

J = Jo exp( / V ■ u dt) 
Jo 

Clearly J = 1, since 0(£, 0) = Thus it is clear that J remains non-zero if 
Given t 



< oo. 

hat g(x) e W 1 ' 1 ^"), there exists Mel, such that 
d 



dx 



9 



< M < oo, 



Vi. 



(30) 



L 1 
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Knowing that -S-g G L 1 and u G L°°, we know that jr-Uj exists and that 



d__ d_ 
dx; 3 dx; 



(31) 



Using this and bounds (26) and (30) , Young's inequality states 



d__ 
dxi J 



< 



d_ 

dxi 



\u 



'3\\L° 



< M \\Uq 



'j I I L° 



(32) 



which leads directly to the bound 

V • udt 



o 



< nM max w . L £. 

3 



(33) 



Thus for finite time, the Jacobian of <j> remains uniquely invertible, with a 



continuously differentiable inverse and thus (27) is a unique C 1 (IR' 1 ) solution 
to (J2TJ , (g , and (|23|) . 

Theorem 4.2 Lei #(x) G W^^IR™) and u (x) e L°°(]R n ), t/ien there exists 
a unique global solution u(x, t) G L°°(]R n ) to i/ie initial value problem (21) 



, (22) , and (23) . 



Proof The proof is the same as for Theorem 4.1 0(^,t) still has a unique 
continuously differentiable inverse, and the solution remains u(x, t) = uo(0 _1 (x, £)), 
but now lacks continuity due to the initial conditions. 



5 Invariants of Motion 

Inviscid and viscous Burgers equation have invariants of motion that are a 
result of their inherent geometric structure. Often these invariants of motion 
can only be realized with additional constraints which impart some physical 
meaning. These constraints are specifically: restricition to a single dimen- 
sion, assuming the velocity is the gradient of a potential function, u = V</>, 
and by the addition of a continuity equation. Both restricting to a single 
dimension and assuming a potential function establishes a curl free veloc- 
ity. Introducing a continuity equation similar to that in the Euler equations 
introduces density. Table [5] shows the various invariants of motions for the 
equations under different constraints. 
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By restricting inviscid Burgers equation to one dimension there are a 
countably infinite number of conserved quantities [47J. When viscous Burgers 
equation is restricted to one dimension, only the wave mass is conserved. 
Similarly, CFB conserves wave mass when restricted to one dimension [35J. 
This result is dependent upon the filter g being an even function, which was 
a criteria set down in section [3] In the special case of the the Helmholtz 
filter, this quantity is also a Hamiltonian structure j35j [36] . 

One dimensional CFB has an additional conserved quantity, total varia- 
tion, that is not found in either inviscid or viscous Burgers equation. CFB 
can be considered a convection equation, with convective velocity u. Since the 
solution is a remapping of the initial condition, the total variation remains 
constant over time. For continuous initial conditions and thus continuous 
solutions, this can be verified directly by noting that T.V.(f) = J \ f x \dx for 
continuous functions. By taking the derivative of ^ and integrating its 
absolute value, one obtains 



0_ 
dt 



J \u x \dx + J sgn(u x ){uu x ) x dx = 0. (34) 
Break the second term into intervals where sgn(u x ) remains constant. u x 



and u are continuous due to Theorem 4.1, so where sgn{u x ) changes sign, 
u x = 0. Thus the second term is zero and J \u x \ dx remains constant over 
time. This was established by Bhat and Fetecau for the Helmholtz filters 



[36], but remains true for any filter satisfying Theorem 4.1 

When assuming a potential function for the velocity, u = V0, in invis- 
cid and viscous Burgers equation wave mass is again conserved. Notice that 
restricition to one dimension trivially implies a potential function for the ve- 
locity. Similarly CFB conserves wave mass with the assumption of a potential 
function. The velocities u and u have different potential functions related 
through the filter g. Again the conservation of wave mass is dependent on g 
being radially symmetric. 

The addition of a continuity equation to inviscid Burgers equation results 
in multiple invariants of motion (17] . It can be easily shown that this addition 
conserves mass, momentum, and kinetic energy. CFB also conserves mass, 
momentum, and kinetic energy when a continuity equation involving the 
filtered velocity is added. 

CFB has another quantity that remains constant over time. As can be 
seen in the existence and uniqueness theorem in section |4| the infinity norm 
of the unfiltered velocity remains constant. Indeed both the maximum and 
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minimum of the initial condition is preserved over time. This was shown for 
the one dimensional case by Holm and Staley [35J. 

It should be noted that almost all of the invariants of motion are found 
using the unfiltered velocity. The filtered velocity, which is to be considered 
the physical quantity, has only wave mass preserved. In order for the filtered 
velocity to accurately model the macroscopic (low wavenumber) behavior, 
energy must be leaving the system in the high wavenumber spectrum. Indeed, 
in numerical simulation in sections [9] and 10 show that energy is leaving the 
filtered velocity at the high wavenumbers. 

Besides wave mass, the filtered velocity does not appear in the invariants 
of motion. It was mentioned in the introduction that the filtered velocity is 
meant as the physically relevant quantity and is to capture the low wavenum- 
ber behavior. Thus it should not capture all of the energy behavior in the 
higher wave modes. The unfiltered velocity is intended to capture such be- 
havior. Therefore the invariants of motion for the unfiltered velocity would 
be more properly compared to those of inviscid Burgers equations, just as 
the invariants of motion for the filtered velocity should be compared to those 
of viscous Burgers equation. It should be noticed that inviscid Burgers equa- 
tion can be thought of as a conservation law for the wave mass. Thus this 
invariant of motion is of primary importance and is noticeably preserved 
in viscous Burgers equations and in both velocities in convectively filtered 
Burgers equation. 



6 Traveling Wave Solutions 

In this section, one dimensional CFB, ^ and ^ , is shown to have a 
traveling wave solution. For the Helmholtz filter, these have already been 
found by Bhat and Fetecau [3^] and Holm and Staley [35]. This section 
generalizes the results to general filters. This traveling wave solution is a 
moving shock that handles an amplitude drop or increase. It satisfies the 
limiting boundary conditions 

lim u(x,t) = u r , lim u(x,t) = itj. (35) 

X— >0O X— > — oo 

Inviscid Burgers is known to have weak solutions 

l +\ \ u t x <ct , , 

= ( Ur x>ct (36) 
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Additional Constraints 


Viscous Burgers 
System Equations 


Invariants of Motion 


One Dimensional 


u t + uu x = vu xx 


/« 


Potential Function 


Ut + u ■ Vu = v A u, 

u = v0 




Additional Constraints 


Inviscid Burgers 
System Equations 


Invariants of Motion 


One Dimensional 


Ut + uu x = 


/ « n ,Vn e Z 


r otennai r unction 


u, + u- Vu = 0, 
u = v0 




Continuity Equation 


u t + u ■ Vu = 0, 
p t + V(pu) = 


f P, f put, jpu-u 


Additional Constraints 


Convectively Filtered Burgers 

System Equations Invariants of Motion 


One Dimensional 


Ut + uu x = 


ju, ju, TV(U), IHloe 


Potential Function 


u t + u • Vu = 0, 

u = v/0 


§ Hit J l^ii | |u| |oo 


Continuity Equation 


u t + u ■ Vu = 0, 
p t + V(pu) = 


/ p, / ptii, / pu • U, ||u||oo 


Table 2: This table lists 


the additional constraints, 


the modified equations, 



and the corresponding invariants of motion for viscous, inviscid, and convec- 
tively filtered Burgers equation. 
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with c = \{u r + ui). The speed at which discontinuities travel is dictated by 
the Rankine-Hugoniot jump condition [T6] . 



It is not difficult to show that (36) is a weak solution to one dimensional 



CFB as well. With Equation (36) as the unfiltered velocity solution, the 
filtered velocity is 

/x— ct 
g a (s)ds + u h (37) 
■oo 

with g a as the filter and c = \{u r + u{). As filter g is an even function, the 
value of the filtered velocity in the middle of the shock is the speed of the 



shock, u(ct,t) = c. This is crucial in verifying that (36) is a weak solution. 

There are some notable properties of the traveling wave solution. It trav- 
els at the same speed as traveling wave solutions to inviscid Burgers. As 
it is desirable for CFB to capture wave speed accurately, this is promising. 
Additionally, as a — > 0, the filter g approaches the Dirac delta function, thus 



(37) converges to (36) , thus the traveling solutions to CFB converges to 
weak solutions of inviscid Burgers equation. 

It is also noted that when u\ < u r that the traveling wave solution will 
converge to a discontinuity, where u(ct~,t) < u(ct + ,t). The Lax entropy 
condition states is that at points of discontinuity u(x~,t) > u(x + ,t). Thus 
the traveling wave solutions to CFB will converge to a weak solution of invis- 
cid Burgers equation which violates the entropy condition. This non-entropic 
behavior can be avoided by only allowing continuous initial conditions. This 
is shown in 



7 Numerics 

For numerical simulations a pseudospectral method with a third order to- 
tal variation diminishing (TVD) Runge-Kutta scheme was utilized. Spatial 
derivatives and filtering were performed in the Fourier space. The equations 
were advanced in time using the optimal third order TVD Runge-Kutta 
scheme presented by Gottlieb and Shu jl9]. The CFL coefficient was chosen 
at c = 0.3, with runs typically conducted at resolutions of 2 12 = 4096 for one 
dimension and 256 x 256 for two dimensions. Aliasing errors were handled 
using the same technique as in Holm and Staley [35J. 

It should be noted that this numerical method captures the behavior in 
both the unfiltered and filtered velocity. The filtered velocity is considered 
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the physically relevant quantity and requires less resolution to capture, thus 
a more efficient numerical method would resolve only the filtered velocity 
using Equation (10) . However, since the priority currently lies in exploring 
the properties of CFB rather than in making efficient simulations, the scheme 
was designed to capture both velocities. 

Our groups has extensive experience in using such techniques with LAE-a 
and LANS-a simulations where the accuracy and stability of the approach 
has been established [5]. Of particular interest, we have shown that the H 1 
norm of the energy in LAE-a computations are numerically conserved in 
longtime calculation [3]. 



8 Shock Behavior 
8.1 Unfiltered Velocity 

In section |4j it was proven that the unfiltered velocity has a continuous 
derivative for all time. This was shown by proving that the Jacobian of the 



material map remains nonzero. In Equation 29 the Jacobian is shown to 
have an exponential structure. As it happens the Jacobian remains nonzero, 
but approaches zero rapidly. From a method of characteristics point of view, 
the characteristics are growing continuously closer, but never intersect. Thus 
the unfiltered velocity will form shocks of continually smaller thickness. In 
numerical simulations, this process will inevitably drop below the finite res- 
olution, thus from a numerical perspective, the unfiltered velocity becomes 
discontinuous. The filtered velocity, however, appears continuous in the nu- 
merical simulations. The shocks in the filtered velocity will not drop below 
the resolution of the numerical run provided a large enough a has been cho- 
sen. For the average run, a is chosen to be at least a magnitude greater than 
Ax. This puts approximately 15 points across the width of the shock, which 
seems a reasonable number to resolve its features. 



8.2 Shock Profile 

In ID viscous Burgers equation, a wave will propagate to the right, with 
the higher velocities overcoming the lower velocities, creating steep gradi- 
ents. This steepening continues until a balance with dissipation is reached. 
Figure [2] shows the evolution of viscous Burgers equation with a Gaussian 
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distribution as an initial condition. 




In figure [3j the evolution of CFB can be seen for two filters. The behavior 
is qualitatively the same as viscous Burgers equation in reference to wave 
propagation and shock formation. Every filter that has been numerically 
simulated has shown similar behavior. 

Numerical simulations were also conducted in two dimensions. It is more 
difficult showing shock formation in 2D, but figures [4] and [5] show the evo- 
lution of a pulse under the Helmholtz and Gaussian filter. In both runs, 
the pulse moves up and to the right, becoming steeper, but never forming 
a discontinuity. Once the shock has formed, the pulse begins to decrease in 
amplitude. 
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(a) 



(b) 



Figure 3: Evolution of CFB using different filters from t — 0, . . . , 10. Only 
the filtered velocity is shown. For all filters a = 0.05. (a) Helmholtz filter 
g(x) = \exp(—\x\). (b) Gaussian filter, g(x) = n~ 1 / 2 exp(—x 2 ). 





(a) 



(b) 



(c) 






(d) 



(e) 



(f) 



Figure 4: Contour plots showing the evolution of 2D CFB using a Helmholtz 
filter with a Gaussian pulse as the initial condition, a = 0.08 with 128 x 
128 grid points. 
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(b) 






(d) (e) (f) 

Figure 5: Contour plots showing the evolution of 2D CFB using a Gaussian 
filter with a Gaussian pulse as the initial condition, a = 0.08 with 128 x 
128 grid points. 



8.3 Shock Thickness 

One of the characteristics of viscous Burgers equation is that as v becomes 
smaller, the shocks formed by the solution become thinner and steeper. 
Heuristically, smaller z/s allow more energy into higher wavemodes, which 
causes a steeper gradient. 

Numerical simulations were run for various values of a and various filters 



in Equations (11) and (12) , to examine the effects upon shock thickness. 
Smaller a correlates with less dampening in the high wavemodes, thus allow- 
ing steeper shocks much like the effect of the viscous term. Figure |6^i shows 
shocks for the Helmholtz filter for different values of a. As a decreases the 
shocks get thinner. Similar results hold true for all filters and simulations 
made in two dimensions. 

Shock thickness can be examined analytically by looking at the traveling 
wave solution. Here the thickness of the shock is defined to be the length 
over which 90% of the amplitude change takes place, centered at the center. 
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X 



(a) (b) 

Figure 6: The thickness of the shocks formed in CFB vary depending upon 
the value of a. (a) As a decreases the thickness of the shock decreases. 

a = 0.08= , a = 0.05= , a = 0.02= (b) The thickness 

of the traveling shock decreases linearly with a. Helmholtz filter, g(x) = 
7,exp(—\x\)= , Gaussian filter, g(x) = ix~ l l 2 exp(—x 2 )=- - - -, Hat fil- 
ter, g(x) = {x — 1 for x G (— 1, 0), 1 — x for x G (0, 1), otherwise}= — ■ — , 
Tophat filter, g(x) = {1 for x G [— |, |], otherwise} = — • • — . 



From section |6j the traveling wave solution is 

/x—ct 
g a (s)ds + ui. (38) 
■oo 

The thickness of the shock will then be 2ab, where b is the value where 

f b f ab 1 x 

/ g{x) dx = / -g{~) dx = 0.9. (39) 

J-b J-ab a " 

This length is independent of u r and U[. As such, the thickness of the shock 
varies linearly on the parameter a. Figure |6Jd shows shock thickness versus 
a for different filters. 



9 Spectral Energy 

Analytically CFB has been examined on an infinite domain. This is of course 
impossible numerically, so in numerical experiments, simulations were per- 
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formed on the domains [0, 27r] and [0, 2tt] x [0, 2tt] for one and two dimensions, 
with periodic boundary conditions. Since simulations were performed with 
a pseudospectral method, obtaining the spectral energy decompositions was 
easily done. 

The examination of the spectral energy decompositions begins by stating 
a special case of Sobelev Embedding Theorem that is found in Hunter and 
Nachtegaele (50]. There it is stated that for a function / : T n — > C defined 
by 

/(x) = ****** ( 4 °) 

meZ™ 

|m| 2p |a m | 2 < oo (41) 

meZ n 



that if 



for some p > § then / is continuous. Furthermore if condition (41) holds 
true for p > j + | then / has j continuous derivatives. 
If spectral energy E{k) is defined as 

E(k) = w 2 , ( 42 ) 

|m|=fc 



then (41) can then be rewritten as 

oo 

k 2p E{k) < oo. (43) 

A:=0 

Thus it can be seen that the rate at which E(k) decays as k — * oo, can 
determine the smoothness of the equation. In one dimension, if E(k) decays 
faster than p, continuity is guaranteed. In a logarithmic plot, this correlates 
with an energy cascade slope of less than -2. In n dimensions, if E(k) decays 
faster than (energy cascade slope less than -{n + 1)) continuity is guar- 
anteed. Existence of continuous derivatives can be guaranteed in a similar 
fashion. 

Inviscid Burgers and viscous Burgers (in its inertial frame) have been 
shown to have an energy cascade slope of -2 during shock formation jHJ 25] . 
Numerical simulations suggest that CFB also has an energy cascade slope 
of -2 for the unfiltered velocity, independent of filter and dimension. The 
filtered velocity's energy cascade slope is then highly dependent upon the 



21 



filter employed. In figure [7j the spectral energy decompositions for simula- 
tions with the Helmholtz and Gaussian filter in one and two dimensions are 
shown. The spectral energy decompositions are shown for a time well after 
the shocks are fully developed. The -2 energy cascade slope is clearly seen, as 
well as the effect of the filters on the filtered velocity's energy cascade slope. 
The Helmholtz filtered velocity displays a -6 energy cascade slope for high 
wave numbers, guaranteeing a continuous derivative in one and two dimen- 
sions. The Gaussian filtered velocity decays faster than any polynomial for 
high wave numbers, thus guaranteeing infinite continuous derivatives. This 
is similar to the solutions of viscous Burgers equation. The viscosity term 
causes an exponential drop in the energy spectrum at high wave numbers, 
also guaranteeing infinite continuous derivatives. 

In section [4], the unfiltered velocity and thus the filtered velocity were 
proven to have a continuous derivative for continuously differentiable initial 



conditions. In section |8.1[ it was suggested that despite this, the unfiltered 
velocity can form shocks that will grow narrower than any finite resolution. 
Thus for numerical purposes the unfiltered velocity can be considered discon- 
tinuous. The filtered velocity will appear continuous to numerical simulations 
if the energy cascade slope satisfies conditions dictated by the special case of 
the Sobolev Embedding Theorem mentioned. 



10 Energy Norms 

The kinetic energy for Burgers equation can be defined as J • u. In CFB, 
there are two different velocities presenting three different kinetic energies: 
J ■ u, J • u, and J |u • u. Analytical comparisons between the energy 
decay rates of viscous Burgers equation and CFB have proven fruitless with 
the exception of the one dimensional case using the Helmholtz filter. In one 
dimension, the nonlinear term is more easily handled, and the Helmholtz 



filter (14) provides a convenient inverse, allowing greater manipulation 
capabilities. 

For one dimensional viscous Burgers equation, kinetic energy can be de- 
fined as 

m = r u^tf^ (44) 
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Figure 7: Spectral energy decompositions for ID and 2D CFB using 
Helmholtz and Gaussian filters . All simulations were performed with 
a = 0.08. (a) and (b) are one dimensional simulations at t = 3 with ini- 
tial conditions u(x) = sin(x) for Helmholtz and Gaussian filters respectively. 
In figure (a) one can see that the energy cascade slope drops dramatically 
after wavelength -. This occurs similarly in figure (b), with the exception 
that the spectral energy decreases exponentially due to the Gaussian filter, 
(c) and (d) are two dimensional simulations at t — 1 with Gaussian pulses as 
the initial conditions, again with Helmholtz and Gaussian filters respectively. 
The spectral energy decompositions show similar characteristics as the one 
dimensional simulations. The spectral energy of u , u - - - -, and a 
reference slope of -2 — • — are shown. 
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The decay rate of the energy is easily calculated to be 

d fu 2 



dx = —v \ u 2 x dx. (45) 



dt 

For one dimensional CFB using the Helmholtz filter, the energy decay rates 
are 



/2 
ycb = a 2 J ( ^{u x + u x )dx (46) 

( ' ^ " U dx = a 2 J(u x ) 3 dx (47) 



d f u 2 
dt 



dt J 2 
d_ fu 2 
dt 



J —dx = —3a 2 J u x u xx {I — a 2 dx 2 ) 1 {u)dx. (48) 



Here a similarity in structure can be seen, especially when comparing (45 ) 



with (47) . The integrands of both are the first derivatives of velocities, 
implying that the majority of the energy is lost through steep gradients, 
a common concept. The similarity is interesting considering that in one 
dimensional viscous Burgers equation, the energy is lost through the viscous 
term, while the energy is lost in the nonlinear term for one dimensional CFB. 

Numerical simulations were conducted for both viscous Burgers equation 
and CFB in one and two dimensions. Figure [8] shows the energy decay rates 
for viscous Burgers equation. Figure [9] shows the energy decay rates for CFB 
using the Helmholtz and Gaussian filters. In figure [9j it can be seen that 
the different kinetic energies tend toward each other with primary difference 
occurring in amplitude. Comparing figures [8] and [9j it is seen that the energy 
decay rates behave quantitatively similar between CFB and viscous Burgers 
equation. Energy decay is minimal until the formation of a shock and then 
the energy decreases rapidly upon formation of the shock. 



10.1 Energy Can Increase 

It is also important to notice that Equations (46) , (47) , and (48) are not 
sign definite as is Equation (44) . This shows while the energy for viscous 
Burgers equation must always decrease, CFB can experience an increase in 
energy. Specifically, it can be seen in Equation (47) , that for steep decreasing 
gradients energy will be lost, but for steep increasing gradients energy will 
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(b) 



Figure 8: Energy decay for ID and 2D viscous Burgers equation, (a) The 
energy decay for ID viscous Burgers equation. Initial conditions u(x, 0) = 
exp(— 3x 2 ), with v = 0.001 and 1024 gridpoints. (b) The energy decay for 
2D viscous Burgers equation. Initial conditions ui(x,y,0) = exp(—30x 2 — 
30y 2 ) and u 2 (x,y,0) = exp(-30x 2 - 30y 2 ), with v = 0.001 and 128 x 128 
gridpoints. J u • u — • — . 



be gained. However, due to the nature of the equation increasing gradients 
will decrease in steepness, while decreasing gradients will increase or remain 
steep. This can be seen in the evolution of a Gaussian pulse, as in figure [3] 
Thus any increase in energy should be brief. 

One would also suspect that as a — > that the brief energy increase would 
dissappear. If solutions to CFB are, in fact, limiting to solutions of inviscid 
Burgers equation, then it follows that the energy rates would approach those 
of inviscid Burgers equation, as well. 

Figure 10 shows two numerical simulations where the initial conditions 
were chosen to generate an increase in energy. In both cases the energy will 
briefly rise before beginning to decay as expected. Figure 11 shows three 
simulations with decreasing a. As was expected, the brief increase in energy 
becomes less dramatic as a decreases. 
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Figure 9: Energy decay for ID and 2D CFB. j u • u — • — , J u • u - - - -, 
j u • u — . Figures (a) and (b) show the energy decay for ID CFB with 
the Helmholtz and Gaussian filter respectively. Initial conditions u(x, 0) = 
exp(—3x 2 ), with a = 0.05 and 1024 gridpoints. Figures (c) and (d) show 
the energy decay for 2D CFB with the Helmholtz and Gaussian filter. Initial 
conditions ui(x,y, 0) = exp(— 30x 2 — 30y 2 ) and U2(x,y, 0) = exp(— 30x 2 — 
30y 2 ), with a = 0.05 and 128 x 128 gridpoints. Decay rates are similar to 
those seen in figure |8j 



11 Conclusion 

By passing the convective velocity through a low pass filter, the cascade of 
energy into higher wave modes of inviscid Burgers equation has been altered. 
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(b) 

Figure 10: Here the energies for CFB are shown, j u • u — • — , J u ■ u 
- - - -, J u • u — . In (a) the initial conditions are uq = C(x — 7r)/(l + 
100(:r — vr) 4 ), with C chosen such that max(iio) = 1, and (b) with a random 
initial condition. In both cases the energies increase initially, but then behave 
with normal energy decaying behavior. 



For a general class of filters, continuous initial conditions lead to continuous 
solutions for all time. Much of the geometric structure seen in invariants of 
motion is preserved in the averaging process. We can see that traveling wave 
solutions propagate at the correct wave speed, and converge to weak solutions 
of inviscid Burgers equation. Through numerical simulations, much of the 
behavior of CFB can be examined. The shock formation and behavior of the 
unfiltered velocity appear similar to that of viscous Burgers equation. The 
shock thickness can be regulated by the parameter a, which is the character- 
istic width of the filter. Spectral energy decompositions give insight into how 
the highwave modes are handled and what filters are needed to guarantee the 
smoothness of the solution. Energy norms also compare favorably to those 
found in viscous Burgers equation. All together, filtering the convective ve- 
locity appears to be a valid technique for high wavemode regularization and 
can hopefully be extended beyond Burgers equation into more general fluid 
dynamics equations. 



a 
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Figure 11: The energy J u ■ u is compared for three different values of a. 

a = 0.08 — • — , a = 0.05 - - - -, and a = 0.02 . The initial conditions 

are uo = C(x — 7r)/(l + 100 (x — 7r) 4 ), with C chosen such that max(-u ) = 1. 
For all three values of a there is a brief increase of energy, but as a becomes 
smaller, this increase in energy becomes less substantial. 
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